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ABSTRACT 

Context. 4U 1538-52, an absorbed high mass X-ray binary with an orbital period of ~3.73 days, 
shows moderate orbital intensity modulations with a low level of counts during the eclipse. Sev¬ 
eral models have been proposed to explain the accretion at different orbital phases by a spherically 
symmetric stellar wind from the companion. 

Aims. The aim of this work is to study both the light curve and orbital phase spectroscopy of 
this source in the long term. Particularly, the folded light curve and the changes of the spectral 
parameters with orbital phase to analyse the stellar wind of QV Nor, the mass donor of this binary 
system. 

Methods. We used all the observations made from the Gas Slit Camera on board MAXI of 4U 
1538-52 covering many orbits continuously. We obtained the good interval times for every orbital 
phase range which were the input to extract our data. We estimated the orbital period of the system 
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and then folded the light curves and we fitted the X-ray spectra with the same model for every 
orbital phase spectrum. We also extracted the averaged spectrum of all the MAXI data available. 
Results. The MAXI spectra in the 2-20 keV energy range were fitted with an absorbed Comp- 
tonization of cool photons on hot electrons. We found a strong orbital dependence of the absorp¬ 
tion column density but neither the fluorescence iron emission line nor low energy excess were 
needed to fit the MAXI spectra. The variation of the spectral parameters over the binary orbit were 
used to examine the mode of accretion onto the neutron star in 4U 1538-52. We deduce a best 
value of M/Voo = 0.65 x 1(T 9 M a yr~ l /(km j - 1 ) for QV Nor. 

Key words. X-rays: binaries - pulsars: individual: 4U 1538-52 


1. Introduction 


4U 1538-52 is a high mass X-ray binary (HMXB) pulsar with a B-type supergia nt companion. 


QV N o r. The orbital pe r iod of the binary system is -3.73 days (IDavison et al 


( 2000 ); 


MuWieriee et al 


(Davison et al 


(1977); 


Clarki 


( 2006 )), and the magnetised neutron star has a spin period of —529 s 


(ll977h :]Becker et al 


tween 4.0 and 7.4 kpc (Crampton et al 


1977)). The distance to the so u rce was estimated to b e be 


1978); 


Ilovaisky et al 


(1979); 


Reynolds et al 


(1992)), and 


we use its latest estimation itAn this work, i.e., 6.4 kpc. Although an orbital low eccentricity was 


estimated to be e -0.08, (Corbet et al. 


orbit for this system (Clark (2000); 


1993 )), a higher va lue of e - 0.17 suggested an ellip tical 


Mukheriee et al 


(2006)). Nevertheless, 


Rawls et al 


(2011) ar¬ 


gued that a circular orbit and an eccentric orbit cannot be distinguished. The X-ray luminosity is 


minimum during the eclipse which lasts -0.6 days (Becker et al. 


(1977)). 


vatories in different energy bands 

such as Tenma 

(Makishima et al. (1987 

)), EXOSAT ( 

Robba et al. 

(1992 

)), Ginea ( 

Bulik et al 

1995 

), RXTE (Clark 

(2000 

)), BeppoSAX 

Robba et al. 

2001 

)), XMM- 

Newton (Rodes-Roca et al. 

2011 

)), INTEGRAL ( 

Temphill et al. 2013 

, and Suzaku 

Hemphill et al. 


2014). The X-ray characteristics derived from these analyses suggest that the mass accretion onto 


the neutron star is consistent with a spherically symmetric stellar wind from the companion star 


(Mukheriee et al 


(2006)). The X-ray spectrum based on a broad band energy range has been de¬ 


scribed by different absorbed power-law relations, modified by a high-energy cutoff including a 
fluorescenc e iron line at 6.4 keV, whenever present, and with cyclotro n resonant scattering fea tures 


at -21 keV (Clark et al 


(1990); [Robba et al 


(2001)) and at -47 keV (Rodes-Roca et al. 


(2009)). A 


soft e xcess component is f ound in the X-ray s pectrum from BeppoSAX and XMM-Newton observa 


tions (Robba et al. 


( 2001 ); 


Rodes-Roca et al. 


(120111) ). Similar to other wind-fed pulsars, the X-ray 


flux changes while the neutron star is moving in its orbit around the circumstellar environment 
of the companion star. Therefore, the variability is associated with the mass-loss rate by stellar 
wind captured by the gravitational field of the neutron star ( M ). It also depends on the distance 
between both components and the stellar wind velocity, mainly. The estimated X-ray luminosity is 
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(2.6-9.1)x l0 36 e rgs' 1 in the 3-100 keV range, assuming an isotropic emission and a distance of 


6.4 kpc (Reynolds et al 


(1992)). Thus, the analysis of the X-ray spectrum from the neutron star at 


different orbital phases provides us the variability of the model parameters we can use to compare 
with accretion models. 

The Monitor of All Sky X-ray Image (MAXI) presents both all sky coverage and moderate 
energy resolution which gives us the possibility to investigate the orbital light-curves and the orbital 
phase resolved spectra of 4U 1538-52. In this paper, we used MAXI data to study the light curves, 
the orbital phase averaged spectrum, and the orbital phase resolved spectra of 4U 1538-52 in the 
2-20 keV energy range. Contrary to previous studies which focus on specific orbital phases at 
specific times, the more than five years of MAXI data, analysed in this work tend to smear out 
short time scale variations and give more weight to long term accretion structures. We describe the 
observations and the data reduction in Sect. [2] present the timing analysis in Sect. [3] the spectral 
analysis in Sect.[4|and in Sect. [3 we summarise our results. All uncertainties are hereafter given at 
the 90% (At 2 =2.71) confidence limit, unless otherwise specified. 


2. Observations 


MAXI is the first astronomical mission on the International Space Station (ISS; 


Matsuoka et al 


(2009)). MAXI attached to the Japanese Experiment Module-Exposed Facility (JEM-EF). It con- 


sists i n two t ypes o f X-ray slit cameras. The main X-ray camera is the Gas Slit Camera (GSC; 


Mihara et al 


(2011 


Slit Camera (SSC; 


)) operating in the 2-20 keV energy range. The second one is the Solid-state 


Tomida et al 


(12011 )1 operati ng in the 0.7-7 keV ene rgy range. The i n-orb it 


Sugizaki et al 


(2011) and 


Tsunemi et al 


( 2010 ), 


performance of GSC and SSC is presented in 
respectively. 

Every hour and a half MAXI covers almost the entire sky per ISS orbit. Therefore, the source is 
observed around 1 ks per day by MAXI. We accumulated the exposure time to extract the spectra 
with enough level of counts. 


3. Timing analysis 


v'Jj I 


We have first obtained the MAXI/GSC 1-orbit light-curve^ of 4U 1538-52 from MJD 55 058 to 
MJD 56 821 (i.e. nearly 5 years) in the 2-20 keV energy band to estimate the orbital period of the 
binary system. Then, we have searched for a_period assuming a sinusoidal signal and the error was 


obtained using the relations derived by 


Leahvl (119871). Our re s ult wa s P or b = 3.7285 + 0.0006 days 


which is consistent with the value given by 


Mukherjee et al 


(2006) (P orb = 3.728382 ± 0.000011 


days). Then we have folded the light curves with the best orbital period to obtain ene rgy re solved 


orbital intensity profiles. The orbital phase reference is taken from 
the phase zero corresponding to mid-eclipse in Fig.Q] 


Mukherjee et al. 


(2006), with 


1 MAXI homepage http://maxi .riken.jp 
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Fig. 1. MAXI orbital phase resolved background subtracted light curve s of 4 U 1538-52 in 2^1 keV, 4-10 
keV, and 10-20 keV energy ranges using orbital parameters reported bv lMukheriee et ah (2006), and folded 
it relative to the mid-eclipse phase. 



The MAXI/GSC data used in our analysis were extracted by the MAXI on-demand process using 
a 1.6 degrees radius circular region centred at the X-ray source position with an annulus background 
region where we excluded other bright sources in the field. In Fig. Q] we show the background 
subtracted light curves for the MAXI/GSC camera, folded at the orbital period of ~3.73 days, in the 
energy ranges 2-4 keV, 4-10 keV, and 10-20 keV. We note that the X-ray flux during the eclipse is 
compatible with no counts in all the energy bands taking the uncertainties into account. 


4. Spectral analysis 

4. 1. Orbital phase averaged spectrum 


We have extracted the orbital phase averaged spectrum of 4U 1538-52 with MAXI/GSC for the 
same observation duration using the MAXI on-demand processing. For spectral analysis we used 
the XSPEC version 12.8.1 (Arnaud (1996)) fitting package, released as a part of XANADU in 
the HEASoft tools. We used FTOOL grppha to group the raw spectra as 2-133 channels by 4, 
134-181 channels by 8, and 182-1199 channels by 12 so that all the spectral bins were Gaussian 
distributed. We tested both phenomenological and physical models commonly applied to accreting 
X-ray pulsars. 

Comptonization models attempt to provide a better physical description of the underlying emis¬ 
sion mechanisms. We fitted the 2-20 keV energ y spe ctrum with several successful models. We 


have chosen CompST (in XSPEC terminology; 


Sunyaev & Titarchuk (1980)) as the representative 


model, i.e., a Comptonization of cool photons on hot electrons modified by an absorbing column 


http://maxi.riken.j p/mxondem 
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Table 1. Fitted parameters for the MAXI spectra in Fig. [2] 


Parameter 


CompST 


Nh (xlO 22 cm 2 ) 4.9+1.4 


kT (keV) 

T 

Norm (xlO -2 ) 
Flux (2-20 keV) 
A 2 /dof 


4 7 + 0.7 
- 0 . 5 , 

17 - 8 !^ 

2.7+0.6 

6 5 +0 ' 3 
u.j_o.2 

40.6/43 


Notes. A CRSF was included in the model (see text). Unabsorbed flux in units of (xlO 10 ) erg cm 2 s 1 . 

along our line of sight. As MAXI/GSC does not have good enough sensitivity to detect the fluores¬ 
cence iron emission line at ~6.4 keV we have not needed to include it in our fits. The absorption 
column Nh fitted with other models were consistent, within uncertainties, and we will use the 
CompST values in the rest of this work. The soft excess at low energies (0.1-1 keV) presents in 4U 
1538-52 and other HMXBs has not been considered here because it lies outside the energy range of 
the MAXI/GSC spectrum. Nevertheless, the well known fundamental cyclotron resonant scattering 
feature (CRSF) at ~21 keV modifies the X-ray continuum above 16 keV (see references in Sect. [I]). 
Therefore, we inc l uded i n our models the CRSF fixing its parameters with the values ob tained by 


Rodes-Roca et al 


(2009). We used the absorption cross sections from 


Verner et al. 


abundances are set to those of 


Wilms et al 


(1996), and the 


(2000). The fitted parameters for the continuum model 


is listed in Table [I] In Fig. [2] we plot our data together with the absorbed CompST best-fit model, 
and residuals of the fit as the difference between the observed flux and model flux divided by the 
uncertainty of the observed flux. 

As the total Galacti c H i column density in the fine of sight of 4U 153 8-52 is in (9.14 - 9.70) x 


10 21 atoms cm 2 range (Dickey & Lockman (1990), 


Kalberla et al. 


(2005)), the neutron star is mov¬ 


ing almost all the orbit in additional absorbing material related to the stellar wind of the companion. 
Therefore, we studied its variability around the orbit extracting phase resolved spectra in the next 
Section. 


4.2. Orbital phase resolved spectra 


Previo u s works in this sen s e have been carried out for only one or two orbits (IMakishima et al 


. 


(1987). Ro bbaetal 


( 2001 ), 


phase (Rodes-Roca et al. 


Mukherjee et al 


(2006), 


Rodes-Roca (2007)) or for a shorter orbital 


(2011)). We have obtained orbital phase resolved spectra of the HMXB 


pulsar 4U 1538- 52 accumulating the 60 s duration scans into s i x ~0,5 day orbital p hase bins 


outside of eclipse (Nakahira et al 


( 2012 ), 


Doroshenko et al. 


(2013), 


Islam & Paul (2014)). 


In this analysis, we have fitted the orbital phase resolved spectra with the same model we used 
in the orbital phase averaged spectrum (see previous Section). Although we also extracted the 
MAXI/SSC spectra, the time coverage of the SSC was too low to study spectral and flux changes. 
We therefore concentrate on the GSC data. For all the spectral models, we kept the parameters free, 
the orbital phase resolved spectra were fitted from energy range 2 keV to 20 keV. The range of 
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Fig. 2. Orbital phase averaged spectrum of 4U 1538-52 in 2-20 keV band and X-ray model continuum 
modified by the fundamental CRSF. Top panel: Data and absorbed Comptonization of cool photons on hot 
electrons. Bottom panel shows the residuals between the spectrum and the model. 


the xl values was 0.7-1.1 for all fits. The data together with the absorbed CompST best-fit model, 
and residuals as the difference between observed flux and model flux divided by the uncertainty of 
the observed flux, are shown in Fig. [3 We obtained six orbital phases to keep a good S/N ratio. In 
Fig. [3 we plot the highest and lowest flux spectra and two sets of residuals. 


The variability of the spectral parameters for the absorbed thermal Comptonization with a 
spherical geometry (CompST in XSPEC) is shown in Figure [3 The Comptoni zatio n parameter, 
y = kT r/fffi,, c 2 ), determines the efficiency of the Comptonization process (Titarchuk (1994), 


Prat et al 


(2008)) and we derived its value from the six orbital phase ranges (see Table0. This indi¬ 


cates an efficient process which corresponds to a moderate accretion rate, being higher in the orbital 
phase ranges just before and after the eclipse than in the mid orbital phase ranges. Moreover, it is 


consistent with similar observations of Vela X-l performed by MAXI (Doroshenko et al 


(2013)). 

The unabsorbed flux along the orbit has no significant evolution, (7.2 - 8.8) x 10 lu erg cm 2 s 1 in 
the fitting energy range (2-20) keV, showing that the accretion rate is quite stable, in the phase bins 
considered. 


4.3. The stellar wind in 4U 1538-52 

One way to study the stellar wind in HMXBs is in terms of the variability of the equivalent hy¬ 
drogen column density throughout the binary orbit. In fact, a simple spherically symmetric stellar 
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Fig. 3. Orbital phase resolved spectra of 4U 1538-52 in 2-20 keV band. Top panel: Selected spectra and 
absorbed CompST best-fit model. Bottom panel shows the residuals for the model. 

Table 2. X-ray luminosity and Comptonization parameter 


Orbital Phase 

Lx CIO 36 ergs L ) 

y 

[0.1 -0.2] 

3.6 

3.8 

[0.2-0.3] 

3.9 

4.9 

[0.3 -0.5] 

4.1 

2.6 

[0.5 - 0.7] 

4.0 

2.2 

[0.7 -0.8] 

3.8 

3.0 

[0.8 - 0.9] 

3.3 

3.7 


wind model may describ e the observed orb i tal dependence of the co l umn density for cer tain range 


of the orbital inclination (Clark et al 


(1994), 


Mukherjee et al 


(2006), 


Rodes et al. 


(2008)). In order 


to obtain an estimation of the stellar mass-loss rate, the equivalent absorption column was derived 
by integrating the wind density along the line of sight to the X-ray source, i.e., combining the sim¬ 
ple spherically symmetric wind model and conservation of mass equations. Then, the hydrogen 
number density at distance r from the donor, rtn(r), can also be estimated through the following 
relationship: 


nu(r) = 


X h M 


(1) 


m H Vo, (1 - R-k/r)! 3 4nr 2 

with mu the hydrogen atom mass, Xu the hydrogen mass fraction, M the mass loss rate, v cxl the 


terminal velocity of the wind in the range 1400-2 800 kms 1 (Abbott 


represents the binary separation (Clark 


for an OB supergiant (Friend & Abbott 


1982), where r = 23.4 R e 


2000), the parameter /3 = 0.8 - 1.2 is the veloci ty gradient 


1986), and R„ = 17.2 II.. (Reynolds et al. 1992 ) the radius 
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Fig. 4. Orbital phase changes in the free spectral parameters fitted in 2-20 keV band. Spectral parameters 
associated with the absorbed Comptonization of cool photons by hot electrons ( compST ) model. Unabsorbed 
flux in units of 1CT 10 ergcirT 2 s _1 . 


of the donor. The resultant Nu arises from the integration of along the line of sight. In Fig. [5] 
we show the equivalent hydrogen column density measured with MAXI as a function of the orbital 
phase. 

The lines represent the model calculations of the absorption column density for different values 
of [5 and M (M Q yr~ l )lvoo (km s~’), in units of 1CT 9 . The red line corresponds to p = 0.8 and 
M/voo = 0.3; the blue line corresponds to p = 1.2 and M/voc = 0.3; the green line corresponds to 
P = 0.8 and M/voo = 3.0; and the turquoise line corresponds to p — 1.2 and M/vc = 3.0. Our best 
least-squares fit is the black line which corresponds to p — 1.2 and M/ vqq = 1.3 . The value of the 


other orbital parameters have been taken from 


Clark (2000) and 


Mukheriee et al 


(2006). As can be 


seen, our data can not discriminate between different values of p. We also note that our estimations 
of M/ Vo,, did not change very much with the eccentricity and orbital inclination. None of the fits 
could describe all the experimental data. Consequently, we looked for the upper and lower values of 
a better fit. From the eclipse egress to the mid-orbit, the models with a lower M /Voo fit the data well 
(red and blue lines), but from the mid-orbit to the eclipse ingress, a higher value of M/v«> (green and 
turquoise lines) is required. This fact suggest that the neutron star is crossing an overdense region of 
the stellar wind and/or some trailing accreting material passing through the line of sight. Moreover, 
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phase 


Fig. 5. Variation of column density versus orbital phase. Red, blue, black, green, and turquoise lines represent 
the wind model for different values of and M/v „. See text for details. 


hydrodynamical simulations show that the interaction between the gravitational field of the neutron 
star, the radiation field of the X-ray source and the stellar wind, may lead to the formation of 


large structures of higher density compared to the unperturbed stellar wind (Blondin et al 


1990). 


These structures are formed mainly in the trail of the neutron star. Consequently, they act as X-ray 
absorbers predominantly in the second half of the orbit. In Vela X-l, an excess of absorption at 
late orbital phases is usually observed. This excess has been successfully modelled by means of the 


above-mentioned hydrodynamical simulations (Manousakis et al 


2012). The fact that Vela X-l 


and 4U 1538-522 are similar systems, since both have a donor of spectral type B0-0.5 I and very 
short orbital periods, imply that these structures are very likely to be present in 4U 1538-522. 


Assuming that the eclipse egress (<p = [0.1 - 0.3]) represents the unperturbed wind, the best 
value of M/\’oc would be in the range (0.3 - 1.0) x 10 9 M Q yr l /(kms '). The median value of 
this range, namely 0.65 x 10 9 M Q yr~ l /(km s 1 ), is in excellent a greem ent with that obtained by 


Clark et al 


(1994). For the terminal velocities reported by 


Abbott (1982) this would translate into 


a mass loss of (0.4 - 2.8) x 10 6 M Q yr 1 for QV Nor. 


Nieuwenhuiizen & de Jauerl ([X990) investigated the dependence of M on the stellar fundamen¬ 


tal parameters mass M, radius R, and luminosity L for a large sample of stars. They presented a 
simple parametrization that gave a good description of observed mass-loss rates over the whole 
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Hertzsprung-Russell diagram. We noted that their formula has different numbers in the abstract 
than in the text which we think it is the correct one. Therefore, we used the following expression: 


M = 9.55 x 10” 


1.24 


M 

Mo 


0.16 


R 

R~e 


0.81 


M 0 yr 


( 2 ) 


The astrophysical data for QV Nor are log(L/L G ) = 5.21 + 0.13, R/Rq = 17.2 ± 1.0, and 


M/M q =20 + 4 (Re ynolds et al. 


1 992 ). By substituting these values in the Eq. [2] the mass-loss 
rate was estimated as M = 4.5 x 10 7 A7 0 , yr 1 , in agreement with the smaller value obtained by 
considering a spherically symmetric wind. 

Our values agree also with the AbbotTs correlation between mass- loss rate and luminosity, 


Clark et al 


(1994), using a Ginga pointed 


M = (0.03 - 1.0) x 10~ 6 M e yr _I (Abbott 1982). Finally, I 
observation, obtained a range M = f0.9 ~ l.9) x 10 6 M 0 yr -1 , for the eclipse egress data, compatible 
with our estimation. Therefore, combining all the previous estimates we can deduce a mass loss 
range for 4U 1538-52ofM = (1,2±0.7)x 1 0 -6 M 0 yr 1 . These values would point to terminal wind 
velocities larger than 2 000 km s _1 , closer to the upper limit (2 800 km s _1 ) of the range derived by 


Abbott (1982). 


5. Conclusions 


From the analysis of the MAXI/GSC light curve, we have estimated the orbital period of the bi- 
na ry sy stem, P or b = 3.7285 ± 0.0006 days, being in agreement with the best value derived by 


Mukheriee et al 


(2006). From the unabsorbed flux, the X-ray luminosity in the 2-20 keV energy 


band was found. Lx = (3.8 ±0.5) x 10 36 erg s _1 . We have investigated the long term orbital variation 
of spectral parameters doing an orbital phase resolved spectroscopy of 4U 1538-52 with MAXI , 
nearly five years of data. A higher value of M/ v cx , is clearly needed to fit the absorption column 
from orbital phases 0.4 into eclipse ingress. Using the Nh values up to orbital phase 0.3, which 
we consider the unperturbed wind, we measured the orbital dependence of the column density 
in the system and of the intrinsic source X-ray flux, and found that it was consistent on aver¬ 
age with a spherically symmetric stellar wind from the optical counterpart of the binary system. 
We have estimated the mass-loss rate range of the early-type B supergiant star QV Nor to be 
M = (1.2 ± 0.7) x 1(T 6 M 0 yr -1 . 


The slight asymmetric accretion column distribution between eclipse-ingress and eclipse-egress 
seen by MAXI was also found by RXTE and BeppoSAX ( Mukheriee et al.lLoOb ). suggesting that a 


trailing wake of material around the neutron star seems to be a permanent structure in this source. 
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